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Abstract: The maximum surplus before ruin is an important index of assets in insurance institutions. 
Considering the important impact of random error factors on the nature of sample paths 
of the surplus process, which essentially increases difficulties in research, we investigate 
the distribution of the maximum surplus in generalized Erlang(n) risk model perturbed by 
diffusion in this paper. We derive a homogeneous integro-differential equation with certain 
boundary conditions, describing the maximum surplus. Particularly, we can deduce explicit 
results as long as the individual claim size is rationally distributed. Except for extending 
a number of results of simple generalized Erlang(n) risk model successfully, our arguments 
are more practical and the results are more delicate. 
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1 Introduction 


Most literature on ruin theory has concentrated on classical risk models, in which claims 
occurr as a Poisson process. Andersen!!! assumed the claim process to be occurred as a more 
general renewal process and derived an integral equation for the corresponding ruin probability. 
Ever since then, random walks and queuing theory have been used to describe a more general 
framework of risk models, under which we can deduce explicit results for waiting times or claim 
severities with the Erlang and phase-type distribution. As the Sparre Anderson risk process can 
reflect the claim occurrences more properly than classical risk models, it is of great practical 
value and considerable importance to study and extend this model. 

A lot of references (see [2-7] and the references therein) made a thorough study of the Erlang 
risk model. Further extension to results concerning generalized Erlang(n) claim waiting times 
have been done by [8-11]. A number of valuable results for survival probability, the surplus 
immediately prior to ruin, the deficit at ruin and the expected discounted penalty function have 
been obtained in above works. 

To make the model more practical, we introduce uncertain factors into risk models. While 
the uncertain factor changed the nature of sample paths of the surplus process, which caused 
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technique difficulties. Hence, the research of risk models perturbed by diffusion are specially 
paid attention in risk theory. Gerber first introduced a diffusion process to express additional 
uncertainty of premium income or aggregate claims. The classical risk model perturbed by 
diffusion has been studied by many scholars, such as [12-16] and so on. 

Li and Garrido!7! considered the expected discounted penalty function for the generalized 
Erlang(n) risk model perturbed by diffusion. In this paper, with the same model, we mainly 
study the distribution of the maximum surplus before ruin, which is an important indicator 
of the assets in insurance institutions. For an insurance company, investigating this maximum 
can not only grasp its ability of withstanding bankruptcy, but also provide an important basis 
for carrying out other businesses, such as the dividend and the investment. 

Consider a time-continuous Sparre Andersen surplus process perturbed by diffusion 


Ni) 
U(t)=u+ct- Y Xi+oW(t), t20, (1) 


i=l 


where u > 0 is the initial reserve and c > 0 is the rate of premium income; X;8 represent 
the claims; the counting process {N(t) : t > 0} denotes the number of claims up to time t; 
{W(t) : t > 0} is a standard Brownian motion that is independent of the aggregate claim 


process 
N(t) 


S(t):= So Xi 


and g > 0 is the dispersion parameter. 

X;s are independent and identically distributed (i.i.d.) nonnegative random variables with 
common distribution P and the density function p. Let py = E(X k) be the k-th moment of 
Xı and 


als) = f ” 7 p(a)dz 


the Laplace transform of the density function p. The counting process {N (t) : t > 0} is defined 
as N(t) = max{k : T} + Tz +--+ +Tk < t}, where the claim waiting times {T;}S2, are iid. 
random variables with a common generalized Erlang(n) distribution, i.e. T; is distributed as 
the sum of n independent and exponentially distributed r.v. 8 


Sn = Vi + Vo +--+ Vn, (2) 


where the V; may have different exponential parameters A; > 0. In the above model, it is 
assumed that {T;}22, and {X;}$2, are independent, and cE(T;) > E(X;), providing a positive 
safety loading factor. 

For u > 0, define T = inf{t > 0 : U(t) < 0}, (co, otherwise) to be the ruin time and 
y(u) = P(T < œ | U(0) = u) to be the ultimate ruin probability. 


2 The maximum surplus before ruin 


In this section, we consider the distribution of the maximum surplus before ruin. For 
0< u < b, let E(u, b) = P(suppcycr U(t) < b, T < œ | U (0) = u) denote the probability that 
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ruin occurs from the initial reserve u without the surplus process reaching level b before ruin. 
Obviously, €(u, b) = 0 for b < u. 

In the following, we can derive that €(u, b) satisfies an integro-differential equation. Let 2 
denote the differentiation operator with respect to u and T the identity operator. 

Theorem 2.1 For 0 < u < b, the probability (u,b) satisfies the following integro- 
differential equation 


ME- - Lge) leu) = [eu zbode + Pew) (3) 


with boundary conditions &(0,b) = 1 and &(b,b) = 0, where P(u) = 1— P(u). 
Proof Fix the number j = 0,1,---,n—1 in the sum S; = Vi + V2 +---4 Vj in (2) with 
So = 0, and define, for 0 < u < b, 


€5,;(u,b) = P( sup U(s) <b, T < co, U(T) < 0| 3; = t, U(t) = u), 


a, j(U, b) = P( sup U(s) <b, T < œ, U(T) =0|S; =t, U(t) = u), 


t<s<T 


E(u, b) = P( sup U(s) <b, T < œ| S; =t, U(t) =u) 
t<s<T 


li 


Es, (u,b) + £a, 5(u, b). 


Similarly to the proof of Theorem 1 in [17], we have that 


((2-£¢- LS emn= [eu-2ve@e+Pw, A 
j=l 

[D (2- £2 - i iaeo = f" tlu- zioaz () 
j=l 3 


Adding (4) to (5), we obtain the desired integro-differential equation (3). The boundary con- 
ditions are obvious: since if U(0) = 0, ruin occurs immediately, and then €(0,b) = 1; similarly 
€;(0,6) = 1, j = 1,2,---,m —1; moreover €;(b,6) = 0, j = 0,1,---,n—1. Then we complete 
the proof. 

Next, for 0 < u < b, define 7° = inf{t > 0 : U(t) > b | U(0) = u} to be the hitting 
time that the surplus process upcrosses the level b, and x(u,b) = P(T > 7? | U(0) = u) the 
probability of the surplus process attaining a given level b without ruin occurring. Noting that 
x(u, 6) =1— E(u, b). 

Corollary 2.2 For 0 < u < b, the probability x(u, b) satisfies the following homogeneous 
integro-differential equation 


I Ge rj ðu 7 x pa) b) = A x(u — 2, b)p(x)dx (6) 


with (0, b) = 0 and x(b,b) = 1. 
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Proof It is easy to obtain the equation (6) by Theorem 2.1 and x(u, b) = 1 — E(u, b). 
Remark 2.3 Similarly, we define x;(u,b) = P(T > 7° | S; =t,U(t) = u), j = 0,1,- ,n— 
1, and then have, for m = 1,2,--- ,n—1, 


wnt) [E-E -i Ze 


jal 


with boundary conditions 
Xm(0, b) = 0, Xm(b, b) = 1, m=0,1,:--,n-1. (7) 
3 Main results 


To solve the integro-differential equation (6) with boundary conditions (7), we consider the 
solution of the following homogeneous integro-differential equation 


B(D)v(u) = A v(u — x)p(x)dz, u0, (8) 


where D denotes the differentiation operator and 


B(D) = Il (z = eee Z D?) = ym" 


j=l 


is a 2n-th order linear differentiation operator. 
By the theory of integro-differential equations!”!, the general solution of the equation (8) is 
of the form 


2n 
v(u) = So n(b)vi(u), u20, (9) 
i=1 
where for i = 1,2, --- ,2n, v;(u) are 2n linearly independent particular solutions of (8) and n; (b) 


are any real numbers. Therefore the solution of (6) with boundary conditions (7) is 


x(u, b) = Fnul), O<usd, (10) 


where m (b), n2(b), +- - ,N2n(b) are determined by the following system of linear equations, for 
m = 0,1,--- ,n—1, (by convention I: =1) 


[Ü E-r- o) [Eno vi(u)]| oO (11) 


j=l 
[i (z z iP <p) [S no] = =1. (12) 
j=l 3 7 i=1 


According to above analysis, we know that the solution of (6) heavily depends on the 
particular solution v;(u), i = 1,2,--- ,2n of (8). For simplicity, let w(u) denote a particular 
solution of (8). In the following, we make a detailed discussion about w(u). 
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3.1 Analysis of the function w(u) 
At first, we do necessary preparation work. By the equation (12) of [17], we know that 


B(s) = p(s), s€C, (13) 


which is the generalized Lundberg equation. Let p1,p2,°-* ,Pn-1> with R(p;) > 0 and pn = 0 
be the n roots of (13). Define the operator T, of an integrable real function f as 


T, f(z) = [co toa rel, Rr)>0, 220. 


x 
This operator has many nice properties, which can be found in [6]. 
Next, replacing v(u) with w(u) in the equation (8) and then taking Laplace transform on 
both sides of (8), we have 


7 d(s) 
w(s)= xr, SEC, 14 
(s) B(s) — p(s) (14) 
where 
2n-1 2n 2n—1 on 
d(s)= Sof Y Bwt) = E w™O) YI Best™ 
j=0 k=j+1 m=0 k=m+1 


is a polynomial of degree 2n — 1. If the p1, p2,°-- , Pn are distinct, by the Lagrange interpolation 


dis) = $ alp,)| Il see] 
j=l 


k=l, k£j Pj T Pk 


and by the theory of divided differences 


B(s) — p(s) 


it 


[He EF p3) { Bl pi, p2: we Pn: S] — pl pr, P2; oe, Pn: S|} 


nm 2n n 
a°” T[i- (8 + ai) 
[Ie -0p (ay TaT Tap], 
where B[p1, p2,--- , Pn, 8] and p[p1, p2,°-* , Pn, 8] are nth divided differences of B and p with 
respect to p1, p2,*** ,Pn,8, respectively. The last step holds due to the property of divided 
difference!l: since B[p1,2,°** , Pn; $] is a polynomial of degree n and the coefficient of s” is 
equal to that of s?” in B(s), which is 


on 
a 
2 (Ii= ri) 
B| p1, P2,°** , Pn, S] can be factored as 


n 
ao” TT, (5 + ai) 
nr ? 
2°(TTia1 4) 
where a1, @2,--: ,@n come in pairs of conjugate complex numbers; there is a close connection 
between the operator T, and divided differences, i.e. 


Bi pr, p2,°°° » Pn; 8] = (—1)" sec, (15) 


TsT pnt pn—1 ey Tp, p(0) = (—1)"B[p1, p2, “oy Pns s]. 
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Consequently, (14) can be rewritten as 


z; n 
a Zi- =1 eee kj rain anii) 


w(s) = 
Ae AD (TT 1 sey a )TsTon Ton-1 -+ To p(0) 


, (16) 


inverting (16), we have that 


w(u) = J w(u — Malady + Go + Seem i (17) 


j=1 
where g(y) = h x lly) = hy *---* hn *l(y) with hi(y) = e7%Y for i = 1,2,--- ,n, x denoting the 
convolution and 


uy) = Wer) p, 


Pn-1°" 


` Tp, p(y) g 


Furthermore, the coefficients are as follows 


2 Ai)d(0) 
gen, ei OTE, ai) 
2 CHi AiJd(p:) 
oT Ten, kei (Pk = PTT (a + pi)] 


"= Mea 6) i Il 


l=1, l4j 


ĉo = — 


i=1,2, -n= 1, 


1 z = 1 1 
d(pr)| poe 
ay — aj 2 r me er pr+a;’ ei 


The equation (17) is a defective renewal equation (see Remark 3.1 below), and the function 
w(u) determined by (17) is a solution of (8). 

Finally, we find 2n linearly independent solutions vilu), i = 1,2,--- ,2n of (8). Since the 
above-mentioned function w(u) is uniquely determined by the derivatives of w(u) at 0, ie. 
w'*) (0), for k = 0, 1,--- ,2n — 1, we can search for particular solutions vi(u), i = 1,2,--- ,2n of 
(8) by determination of the initial conditions. 

From [17], we derive that the survival probability ¢(u) = 1 — w(u) satisfies the equation (8) 
and ¢'(0) # 0 (see Remark 3.2 below). Since 4(0) = 0, we take v2(u) = ¢(u) as a particular 
solution of (8). The other particular solutions can be found by specifying the initial conditions 


as 
v (0) = I{k =i- 1}, k=0,1,--- ,2n—1, i = 1,3,- ,2n, 


where I{-} denotes the indicator function. It is easy to prove that the functions vilu), i = 
1,2,--- ,2n defined as above are linearly independent. 
y (17), we show that v;(u), i = 1,3,--- , 2n all satisfy the defective renewal equation 


n-1 


uu) = f nu- Hanan + gat 3 Geer + Yon se a 


j=l 
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where g(y) is given as above and the coefficients are as follows 


2” (Ti: 0s) Bi 


io = eS e 
aoa (TIZ i pe) (Ts ai) 
n ; 2n—i 
pa zA à i) Vito Bmti Pk , k=1,2,---,n-1, 
a% prl IIS, jen (Pi — Pr) TTi (ar + px) 
2 A) n ] n 2n—i EA n 1 1 nps 
ij = B. i P: j=1,2, n. 
d an a ae Xd a ial M, Spd eer 


3.2 The special case of w(u) for rational claim size 

In this section, we consider a special case of the function w(u) when the claim size is 
rationally distributed. If the density function p is rationally distributed, then w has a rational 
Laplace transform that can be inverted explicitly by partial fractions as follows. We assume 
that claim size X; is rationally distributed, i.e. 


— Qm-1(S) 


p(s) = Onley’ R(s) € (hx, 00), (18) 


where m € N*, hx := inf{s € R : Ele-**] < co}, Qm is a polynomial of degree m with 
leading coefficient 1, Qm—1 is a polynomial of degree m — 1 or less, and Qm and Qm_ do 
not have any common zeros. Moreover, since f(s) is finite for all s, with R(s) > 0, equation 
Qm(s) = 0 has no roots with positive real parts. 

Substituting (18) into (14) yields 


4(8)Qm(8) 


Ol = BGIOn@) O24 


sé, (19) 


where B(s)Qm(s) — Qm-1(8) is a polynomial of degree 2n + m with the leading coefficient 


2n 


iS YS E 


Then it can be factored as 


B(8)Qn(8) — Qnal) = -D'r v [Te- oJ Th e+], 


i=1 


where 1, p2,*-* , Pn—1, with R(p;) > 0 and pn = 0 and —Ry, —R2, -++ ,—Rn+m, with R(R;) > 
0, are all the roots of the equation B(s)Qm(s) — Qm-1(8) = 0. If Pi; P2;°**;Pn—1 and 
Ri, Ro,-++ ,Rnym are all distinct, by partial fractions, the equation (19) can be rewritten 
as 


als) = 24 a ra ra (20) 
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where 
eee 
o(a p(T R; 
2” (TTi As) d (91) Om (pi) 


Qi = — + , i= 1,2,- n- 1, 
o?r Jir (Ri + eJl Thar, kilp — pi) 


j = 2” (TTi M)d(—R;)Qm(-R3) a ee en 
pj a?r) Jia (Ry + DETE rez (Re — Ri) TE ed a 


Taking inverse Laplace transform on the both sides of (20) yields 


ag = 


n+m 


w(u) =at Fae + D Bye Fi", u>0. (21) 


In the following, we find 2n linearly independent solutions v;(u), i = 1,2,--- ,2n of (8) as 
given in Section 3.1. By (19), we have that 


di(s)Qm(S) 


6:(s) = B(s)Om(s) — Omals) =O, se, (22) 


where 
2n—i 


2n 
di(s) = Sa = Se sl. 
k=i t=0 
Making similar discussion as w(u), we obtain 


2” (T à) a. 23 
o [Tan (oy — AG (s +R) 


If p1, P2,°** )Pn—1 and Rj, Re,--- ,Rn4m are all distinct, by partial fractions, we obtain 


;(s) = 


n— n+m 
Bij 


z Qi, 0 Qi, k ij 
; = —— 24 
ĉ:(s) 8 De er s+ R;’ (24) 


where 


2"( TT 1 ri) Bi Qm(0) 


See aa Tye phe Ie). 

ane 2” ( TTi-1 i) di (On) Qm (Pr) karadi 
: om TIPET (Ry + pe) Trea, 1¢n(Pt — Pr) 

Bij = 2 (Lier eal K5) j=1,2,-+>,n+m. 


o| Tia (Ry + PDE E ry (Be -R 


Inverting (24), we yields that 


n+m 


nu) = 040+ ane t2 ofr Rju 4=1,3,---,2n, u20. 
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Remark 3.1 According to the Lemma 1 of [17] and by simple calculations, we obtain that 
the function g(y) is a defective density function 


(Ts pT hes ai) 
which means that (17) is a defective renewal equation. 


Remark 3.2 From the conclusion of Theorem 1 in [17], the survival probability olu) 
satisfies 


i gty)dy = 1 — 


da) = [TI (2- En- FP) aw, m=1,2,---,n—1. 


=} 


ee. 


Let m = 1, we have 


alu) = (I~, D- <D?) u u). 


Taking Laplace transform on both sides of above equation yields 
, = 21 ^ 2 2c 2A1 
P0) = Pails) + (+ Ss- Ty) A). 
When s is large enough, the right side of above equation is greater than zero, i.e. (0) > 0. 
4 The result of (u,b) for n = 2 and à; = 42 =A 


In the following, we consider the Erlang(2) risk model that is perturbed by diffusion, and 
obtain the expression of x(u, b) in terms of the survival probability (u). 
Xing and Zhang!!®) gave that 


Oua(s — p(s +p+ 3 


A) = EBON a 
1 7 AN? O11 
GAD) o2p(a2p + 2c) (26) 


where p > 0 is the unique positive root of Lundberg equation (13), and 


_ 2e — Apa 
0= Nix Apa >0 


is the relative safety loading factor. From Section 2 of [18], it is straightforward to calculate 
the second derivative at zero, i.e. 


BA2 CO u1 2c 
A 0 = a $ i 
o” (0) splotp + 20) 53? (0) 
Using the formula (22) and letting 
2c 
B =p F ee 
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we have 
5 04 p34(s) 
Vv. 8) = OO 
(9) = Debus — per)’ 
which can be inverted as 
USA J jupi f glu =y Pidy, (27) 
0 


where 
An 7 PB 
4761 (p + B) 
From the discussion in Section 3, one knows that 


4 
x(u,b) = $ ni(b)vi(u), 
i=1 


where the coefficients n (b), n2(b), n3(b), n4(b) are determined by the boundary conditions (11) 
and (12). We can solve the values of 7;(b), for i = 1,2,3,4, 


m(b) = n3(b) = 0, 


07(p + B)d(b) + 07 pBlf> o(b — y)e”dy — J? olb — y)e~F¥ dy] 
0? (p + 3)$(b) + [02 p34(b) — 2cd’(b) — o2" (LIS o(b — y)erudy ~ fÈ (b — y)e-Pudy]’ 


n2(b) = 


-[2c¢' (b) + 074"(b)] 
0? A(p + 8)$%(b) + Alo? p8o(b) — 2eg’(b) — 076" (ISa Hlb — y)ervdy — fy o(b — y)e-Pudy] 


Consequently, x(u, b) can be expressed in terms of the survival probability ¢(u). There is the 


na(b) = 


analogous conclusion in Erlang(2) risk model that is given in [7]. As we know, the famous 
Dickson formula holds in the classical risk model, i.e. (u,b) = ga., This means that x(u,b) 
is obtained if we have a knowledge of ¢(u). 

We remark that when claim sizes are exponentially distributed, there is a closed expression 
for ¢(u) (see [17]) and then for x(u, b). 
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EFRA N Erlang(n) MBL AP SRE AE 
TERR, KEE 
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